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PREFACE 


‘This manual presents Mathematical Subroutines which are 
used to assist both the RCA 110A programmer and machine 
operator in the writing and debugging of programs. Each 
subroutine includes a description, a flow chart, and a 
symbolic listing. 


‘The fixed-point routines can be called and executed inde- 
pendently of each other using the SLAP 2 System, or may 
be assembled with the program. 


The floating-point routines may also be called and executed 
by the SLAP 2 System, but care should be taken to call each 
routine needed to run the desired routine. As an example, in 
order to use the Floating Point Common Logarithm both FLG 
‘and DVD would have to be called. 


These routines will be periodically revised and supplemented 
to increase the working tools available to the RCA 110A user. 


vil/ viii 


SECTION 1 
SIN/COS-FIXED POINT 


REVISIONS 
DESCRIPTION. [ DATE | APPROVED 








FIRST MADE FOR 
USED ON 


NEXT ASSY 














CONTRECTENG RADIO CORPORATION OF AMERICA 
NAS 8- 13007 NEW YORK, NY—VAN NUYS PLANT 

DRAWN DATE 

Juv = PL ASE ES SATURN GROUND COMPUTER SYSTEM 


ao a MATHEMATICAL ROUTINE 
ZY OC Dern 1x0 oF SIN/COS- FIXED POINT 


DESIGN ACTIVITY APPD REL] SIZE | CODE IDENT NO. 
2/8 6634 








WEIGHT SHEET 1 0F /2 





DEL 1263-1 (2-65) 


LIST OF MATERIALS OR PARTS LIST 
QTY REQD PART OR % VENDOR ITEM — SEE SOURCE CONTROL OR SPEC CONTROL DWG 
04 ]503 1502 . IDENTIFYING NO. NOMENCLATURE OR DESCRIPTION SPECIFICATION 
2186653 ‘SATURN GROUND COMPUTER SYSTEM 
SPECIFICATION FOR SATURN V 
COMPUTER PROGRAM 
2186634-1 SIN/COS = FIXED POINT 
‘SYMBOLIC CARD DECK 
2186651 SATURN GROUND COMPUTER SYSTEM 
SLAP - 2 SYSTEM 
CDV 86634 CEI DETAIL SPECIFICATION 
MATHEMATICAL ROUTINE 
SIN/COS - FIXED POINT 
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‘SIZE | CODE IDENT NO. 


al49671| e/8 6634 








DEL 1263—2A 


SATURN GROUND COMPUTER SYSTEM 


MATHEMATICAL ROUTINE 


SIN/COS - FIXED POINT 


Prepared Under Contract NAS 8-13007 
For 
National Aeronautics and Space Administration 
Marshall Space Flight Center 
Huntsville, Alabama 
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Mathematical Method 


Program Operation 
Program Assembly and Loading 
Program Control 
Control Program 
Calling Sequence 


Output Format 


Memory Requirements 
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Appendices 
Flow Charts 


Tables 
Program Listing 
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UNLESS OTHERWISE SPECIFIED: DIMENSIONS | SIZE | CODE IDENT NO. 
‘ARE IN INCHES AND INCLUDE THICKNESS 
OF PLATING. TOLERANCES ARE: 


Ka ot 2K = 
FRACTIONS ANGLES = 21/22 
DEL 1263-3 (2/65) 











Introduction 


The purpose of the fixed point Sin/Cos subroutine is to 
calculate the sine or cosine of an argument X where 
-360°€ X <€ 360°. 


This mathematical routine has been prepared to fulfill 
the requirements of Specification Number CDV 86634. 

This drawing, together with the reference documents, de- 
scribes the procedures for preparing program card decks 
and the method of utilizing same. 


The program flow chart shown in Appendix A describes 
the logic used to develop the program. 


Mathematical Method 


The following expression was used to evaluate the Sin X. 


1.5707963 
~ -64596371 
+ 07968968 
~ -00467377 
+ 00015148 


SIN (X +90°) 


Maximum absolute error is less than 10 


RADIO CORPORATION OF AMERICA 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS | SIZE | CODE IDENT NO. R 

ARE IN INCHES AND INCLUDE THICKNESS 

OF PLATING. TOLERANCES ARE: 

Kat 20K 

FRACTIONS ANGLES = 1/22 
DEL 1263-3A (2/65) 



























Program Operation 


Program Assembly and Loading 






This program is assembled by the SLAP 2 System (2186651) 
using the symbolic card deck for this program. The SLAP 2 
Loader uses the object program card deck, a product of the 
SLAP 2 assembly, to load the program into core memory. 
Program Control 
This program is used as a subroutine by a control program. 


Control Program 


Calling Sequence 








Location OP Code Address 





a IDL Argument (see below) 
atl TSP SIN or COS (see below) 
at2 Return 




















Argument: 





The fixed point argument in fraction of a circle 
is scaled at 0 (binary point between bits 23 and 


22) 
SIN - Symbolic entrance in the sine sub-routine 
cos - Symbolic entrance in the Cosine sub-routine 


Output Format 






The result of the calculation will be in the left (L) 
accumulator scaled at 1 (binary point between bits 22 and 
21). 






Error Alam - None 


Memory Requirements 











See Program listing 
Index Registers used - None 
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CODE IDENT NO. 


49671| 2/8 6634 
EC RS CaS Ca 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS 
[ARE IN INCHES AND INCLUDE THICKNESS 
‘OF PLATING, TOLERANCES ARE: 


Mat OK = ae 
FRACTIONS ANGLES = + 1/2° 
DEL 1263-3A (2/65) 












Sin/Cos - Fixed Point Card Deck 


The symbolic source card deck for this program is 2186634-1. 


Appendices 


Flow Chart 


Flow chart for this program is in Appendix A. 


Tables 


Program Listing 


A list of this program showing memory locations (relative 
to zero) octal representation and source card contents 

(ia symbolic language format) for each memory location used 
by the program is in Table 1. 
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AALS cones aap mereve cess 
A thang ToleRANGES ARE 49671| 2/8 663F 
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APPENDIX A 


SATURN GROUND COMPUTER SYSTEM 


SIN/COS - FIXED POINT 


FLOW CHARTS 
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UNLESS OTHERWISE SPECIFIED: DIMENSIONS | SIZE | CODE IDENT NO. 
mene mee 2/86634 


at 


FRACTIONS = * ‘SHEET 


DEL 1263-3A (27/65) 





SET UP EXIT 
ADD 90° TO ARGUMENT 


SHIFT LEFT +1 AND 
COMPLEMENT. 





UNLESS OTHERWISE SPECIFIED: DIMENSIONS | SIZE | CODE IDENT NO. 

See aries wewuoe Tmeuess 

Cregeen A|49671| 2/86 634 
na 


FRACTIONS = + 
DEL 1263-34 (2/65) 





x 2K-1 FOR SIN OF 


ARGUMEN 


UNLESS OTHERWISE SPECIFI IMENSIONS | CODE IDENT NO. | 

ae er arco ae 
rime 49671| 2/8663F 
Kes fa 
nana 

DEL 1263-3A (2/65) 











TABLE I 


SATURN GROUND COMPUTER SYSTEM 


SIN/COS - FIXED POINT 


PROGRAM LISTING 
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(OF PLATING. TOLERANCES A\ 


Kat 
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DEL 1263-30 (2/65) 
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‘SHEET 


2/86 634 


CODE IDENT NO. 
49671 





ANGLES = +1/2° 


DEL 1263-34 (2/65) 





SECTION 2 
SQUARE ROOT-FIXED POINT 


REVISIONS: 


ZONE[LTR DESCRIPTION DATE | APPROVED 
A | SH 5; SECTION 2.0, LINE 3 WAS "WHERE 1> £ [#64 


Feng 
>.25". LINE 7 WAS Yoyp 7 V2 G+ Fy, ) flees oe [Pr thea 


SH 10; BOX 1 OF FLOW CHART, "POSITION" WAS 
“POSITIVE". EQUATION IN BOX 4 WAS Y; = F 
+ 1/4 - 1/2 (F-1/2). EQUATION IN BOX 6 WAS ® 


Yoat 1/2 wy, + F/Y,). 








CLASS II, DCS 642 





[_useo on | 
110A 





FIRST MADE FOR 


NEXT ASSY 














CONTRACT NO RADIO CORPORATION OF AMERICA 
NAS 8-13007 NEW YORK, NY—VAN NUYS PLANT 

DRAWN DATE 

SATURN GROUND COMPUTER SYSTEM 

es cae| MATHEMATICAL ROUTINE 


AO Ze >1nvees} SQUARE ROOT - FIXED POINT 


DESIGN ACTIVITY APPD Ret | SIZE | CODE IDENT NO, 
296711 2/86635 





SHEET 1 OF /# 





DEL 1263-1 (2-65) 


LIST _OF MATERIALS OR PARTS LIST 





QTY REQD PART OR 
IDENTIFYING NO. 


2186653 


% VENDOR ITEM — SEE SOURCE CONTROL OR SPEC CONTROL DWG 





NOMENCLATURE OR DESCRIPTION 
ISATURN GROUND COMPUTER SYSTEM 


SPECIFICATION 





[SPECIFICATION FOR SATURN V 





[COMPUTER PROGRAMS 





2186635-1 


ISQUARE ROOT - FIXED POINT 





[SYMBOLIC CARD DECK 





2186651 


ISATURN GROUND COMPUTER SYSTEM 





ISLAP-2 SYSTEM 





CDV86635 


[CEI DETAIL SPECIFICATION 





[MATHEMATICAL ROUTINE 





ISQUARE ROOT - FIXED POINT 

























































































RADIO CORPORATION OF AMERICA 





DEL 1263—2A 








‘CODE IDENT NO.] 














SATURN GROUND COMPUTER SYSTEM 


MATHEMATICAL ROUTINE 


SQUARE ROOT - FIXED POINT 


Prepared Under Contract NAS 8-13007 
for 
National Aeronautics and Space Administration 
Marshall Space Flight Center 
Huntsville, Alabama 
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UNLESS OTHERWISE SPECIFIED: DIMENSIONS | 
[ARE IN INCHES AND INCLUDE THICKNESS 
OF PLATING. TOLERANCES ARE: 
Mae 20K = 
FRACTIONS ANGLES = 1/2° 
DEL 1263-3A (2/65) 
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INTRODUCTION 


MATHEMATICAL METHOD 


PROGRAM OPERATION 
Program Assembly and Loading 
Program Control 

Control Program 
Calling Sequence 


OUTPUT FORMAT 
MEMORY REQUIREMENTS 
SQUARE ROOT - FIXED POINT - CARD DECK 


APPENDICES 
Flow Charts 


TABLES 
Program Listing 
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uaa ooetes ee teaaes 
Aa IN Wea AD NUDE T agent 218 6635 
Wee peril 
ae ieee yas 
DEL 1263-3A (2/65) 





INTRODUCTION 


The fixed point Square Root routine extracts the square 
root from the 23 most significant bits of an argument in 
the left (L) and right (R) accumulators. The argument is 
initially normalized by the Square Root routine. 


This mathematical routine has been prepared to fulfill the 
requirements of Specification Number CDV 86635. This 
drawing, together with the reference documents, describes 
the procedures for preparing program card decks and the 
method of utilizing same. 


The program flow chart shown in Appendix A describes the 
logic used to develop the program. 


MATHEMATICAL METHOD 
Assume that the argument X is normalized 


Unpack X 


xi/2 2e 1/2 





. £/° where L> £ > .25 


The first fractional approximation, Y,, to el? 46 given 


by £ + 1/4 - 1/2 (£ - 1/2)". ‘This is then followed by 
three Newtonian iterations of the form: 


= 


Yoo. 7 V2, + £/Y,) 


Maximum relative error is less than 107°, 





PROGRAM OPERATION 


Program Assembly and Loading 


This program is assembled by the SLAP2 System (2186651) 
using the symbolic card deck for this program. The SLAP 2 
Loader uses the object card deck, a product of the SLAP 2 
assembly, to load the program into core memory. 


RADIO CORPORATION OF AMERICA 


UNLESS OTHERWISE SPECIFIED: DIMENsions | SIZE | CODE IDENT NO. 
ARE IN INCHES AND INCLUDE THICKNESS 2 /8 6 6 35 
OF PLATING. TOLERANCES ARE: 
ae Ks 
FRACTIONS = ANGLES == 1/22 
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Program Control 


This program is used as a subroutine by a control program. 


Control Program 
Calling Sequence 
LOCATION 


OP CODE ADDRESS. 





a - Argument 


atl 
a+2 


TSP SQR 


Return 


(see below) 


(see below) 


The fixed point argument is in the left (L) 
and right (R) accumulators, scaled at an even 
number, Q. 


Argument: 
SQR Symbolic entrance in the Square Root subroutine. 
OUTPUT FORMAT 


The square root of the argument will be in the L, R 
accumulators scaled at Q/2. 


Error Alarm: If the argument is negative, the square root 
of the absolute value appears in the L, R accumulators 
scaled at Q/2 and the overflow indicator is set. 


MEMORY REQUIREMENTS 


See Program listing 
Index Registers used - none 





SQUARE ROOT - FIXED POINT - CARD DECK 





The symbolic source card deck for this program is 2186635-1. 


RADIO CORPORATION OF AMERICA 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS | SIZE | CODE IDENT NO. 


REIN INGHES AND INCLUDE THCRNESS 
A| 49671 


ane mene poe 2186635 


Mae OK 
FRACTIONS =:t  ANGLES=+1/2° 
DEL 1263-3A (2/65) 





APPENDICES 


Flow Chart 


Flow Chart for this program is in Appendix A. 


TABLES 


Program Listing 


A listing of this program showing memory locations (relative 
to zero), octal representation and source card contents 

(in symbolic language format) for each memory location 

used by the program. 
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UNLESS OTHERWISE SPECIFIED: DIMENSIONS J SIZE | CODE IDENT NO. 
[ARE IN INCHES AND INCLUDE THICKNESS 2 / 8 6 6. 3. Ss 
OF PLATING. TOLERANCES ARE: 
Met Kat 
FRACTIONS ANGLES = 1/2° 
DEL 1263-3 (2/65) 








APPENDIX A 
SATURN GROUND COMPUTER SYSTEM 
SQUARE ROOT - FIXED POINT 


FLOW CHARTS 





RADIO CORPORATION OF AMERICA 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS ‘CODE IDENT NO. 
‘ARE IN INCHES AND INCLUDE THICKNESS Z 18 6 635 
OF PLATING. TOLERANCES ARE: 
Mae 2K 
FRACTIONS =x ANGLES= 1/29 
DEL 1263-3A (2/65) 





CLEAR OVERFLOW INDICATOR 
ZERO OVERELOW SET LOCATIONS 


COMPLEMENT ARGUMENT 
LOAD OVERFLOW SET LOCATIONS 
WITH 20000000 


NORMALIZE ARGUMENT TO 
SATISFY CONDITIONS 1>F > .25 
AND STORE. 





DIVIDE THE NUMBER OF SHIFTS 
TO NORMALIZE BY 2 AND SAVE 





RADIO CORPORATION OF AMERICA, 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS | SIZE | CODE IDENT 
ARE IN INCHED AND WICLUDE THICKNESS 49671 Ea 
met 2K = 


rections == motes=21/2 [SOME aes Cee 


DEL 1263-3A (2/65) 





SHIFT NORMALIZED ARGUMENT 
RIGHT 1 POSITION 


CALCULATE FIRST APPROXIMATION 
Yy = £+ 1/4 - 1/2 (£-1/2) 


CALCULATE: 


You. 7 1/2 cm + £/¥,) 


DECREMENT XR6 BY 1 


IS XR6 = 0? 


RADIO CORPORATION OF AMERICA 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS 
[ARE IN INCHES AND INCLUDE THICKNESS 
OF PLATING. TOLERANCES ARE: 


20K se 
ANGLES = 1/2 





SCALE FINAL ANSWER AT Q/2 


IF ARGUMENT WAS NEGATIVE 
SET OVERFLOW 
LOAD L_& R WITH RESULTS 


RADIO CORPORATION OF AMERICA 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS | SIZE | CODE IDENT NO. 
Grrumnas Tourmncs ac «| A | 29671 2/86635 
mae kee 


mcvons=s mousse [SAE ————S—S—S——i fet 77 


OEL 1263-34 (2/65) 





TABLE I 
SATURN GROUND COMPUTER SYSTEM 
SQUARE ROOT - FIXED POINT 


PROGRAM LISTING 
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OF PLATING. TOLERANCES ARE: 
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FINAL RESCALING 
SET @VERFLOW INDICATOR 
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SECTION 3 
EXPONENTIAL (2, e, 10)-FIXED POINT 


REVISIONS 
ZONE[LTR DESCRIPTION DATE | APPROVED 








A | SH 11, ADDED EXP-3A TO ENTRANCE TO BOX 2, |’““*#4 
ADDED FLOW CHART SH NO. EXP-3. Pres td 
CLASS II, DCS 642 





FIRST MADE FOR 


NEXT ASSY 
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SATURN GROJND COMPUTER SYSTEM 





SLAP - 2 SYSTEM 
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CEI DETAIL SPECIFICATION 





MATHEMATICAL ROUTINE 





EXPONENTIAL (2,e,10)=FIXED POINT 
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for 
National Aeronautics and Space Administration 
Marshall Space Flight Center 
Huntsville, Alabama 
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INTRODUCTION 


The purpose of the fixed point Exponential routine is to calculate 
the exponential to the base 2, e or 10 of an argument, X. 


This mathematical routine has been prepared to fulfill the require- 
ments of Specification Number CDV 86636. This drawing, together 
with the reference documents, describes the procedures for preparing 
program card decks and the method of utilizing same. 


The program flow chart shown in Appendix A describes the logic used 
to develop the program. 


MATHEMATICAL METHOD 


The argument is first converted to its base 2 equivalent using the 
identity: 
ie ot EE, ew os ie 
The result is then shifted to a scaling of 23 to obtain the form 
otf, if 
=z 2 2 


where i and £ are the integral and fractional parts of the argument. 


The following polynoninal is used to compute 2°* 


[ie 


2d = .34657210 


c. 

C, = .06006622 

C, = .00691806 
= .00061973 

C, = 000033177 

Cc, = .000004208 


%: 
| Cc, = 1.0000000 


6 
and i is used to scale the result since 2+ may be accomplished by 
a shift in a binary machine. 


For Xz 0, the relative error is no greater than 10% in magnitude, 
For xe. 0, the absolute error does not exceed 1076 
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PROGRAM OPERATION 


Program Assembly and Loading 


This program is assembled by the SLAP 2 System (2186651) using 
the symbolic card deck for this program. The SLAP 2 Loader uses 
the object card deck, a product of the SLAP 2 assembly, to load 
the program into core memory. 


Program Control 


This program is used as a subroutine by a control program. 


Control Program 
Calling Sequence 


Location OP Code Address 





a - Argument (see below) 
atl TSP Entry Points (see below) 
a+2 Return 


Argument - The fixed point argument is in the left (L) accumulator. 
The scaling of the argument is in the right (R) accumu- 
lator scaled at 23. The approximate maximum values 
are listed below: 


Base X Max 


2 22.99999 
e 15.94238 
10 6.923689 





Entry Points: The program has three entry points. Select one 
of the symbolic entry points listed below depend- 
ing upon the base desired. 


Entry 
EX2 
EXE 
EX10 
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OUTPUT FORMAT 


The result is in the L, R accumulators scaled at 23. 
Error Alarm: If the value of the argument is such that its ex- 
ponential value exceeds the capacity of the L accumulator, the 


overflow indicator is set and the L, R accumulators contain 
erroneous data. 


MEMORY REQUIREMENTS 


See Program listing. 
Index Registers Used - None 


EXPONENTIAL - FIXED POINT - CARD DECK 





The symbolic source card deck for this program is 2186636-l. 
APPENDICES 

Flow Chart 

The flow chart for this program is in Appendix A. 

TABLES 


Program Listing 


A listing of this program showing memory locations (relative to zero), 
octal representation and source card contents (in symbolic language 
format) for each memory location used by the program is in Table 1. 
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SAVE SCALE FACTOR, CONVERSION 
FACTOR AND ARGUMENT BASE 2 
AND_SET_UP_RETURNS 






SAVE SCALE FACTOR, CONVERSION 


FACTOR AND ARGUMENT BASE e 
AND SET_UP_ RETURNS 













SAVE SCALE FACTOR, CONVERSION 
FACTOR AND ARGUMENT BASE 10 
AND SET UP RETURNS 






IS ARGUMENT SCALED AT 
GREATER THAN 21, 4? 


YES 






SCALE ARGUMENT AT 21 
CONVERT ARGUMENT TO BASE 2 AT 23 


AT 23 
‘SIZE | CODE IDENT NO. a 
a\49671| 2/86636 | 


[sce «dT 
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IDENTIFYING NO. 


2186653 


* VENDOR ITEM — SEE SOURCE CONTROL OR SPEC CONTROL DWG. 





NOMENCLATURE OR DESCRIPTION 
SATURN GROUND COMPUTER SYSTEM 


SPECIFICATION 





SPECIFICATION FOR SATURN V 





COMPUTER PROGRAMS. 





2186637-1 


ARCTANGENT ROUTINE ~- FIXED POINT 





SYMBOLIC CARD DECK 





2186651. 


SATURN GROUND COMPUTER SYSTEM 





SIAP-2 SYSTEM 





CDV 8663: 


CEI DETAIL SPECIFICATION 





MATHEMATICAL ROUTINE 








ARCTANGENT ROUTINE - FIXED POINT. 
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Introduction 


The purpose of the fixed point Arctangent routine is to 
calculate the arctangent z. 
This mathematical routine has been prepared to fulfill the 
requirements of Specification Number CDV 86637. This drawing, 
together with the reference documents, describes the 
procedures for preparing program card decks and the method of 
utilizing same. 


The program flow chart shown in Appendix A describes the logic 
used to develop the program. 


hematical Method 


The relative magnitudes of the two arguments are first 
determined. If the absolute value of the contents of A 
is less than the absolute value of the contents of B, then 


us A is found and the arctangent of u computed. If the 

oe @) 
(@) 
is found and the arcotangent of u is computed. The answer 
at this point is in either the first or fourth quadrant. 
The sign of the denominator of the input variable is then 
inspected and, if negative, 180 degrees are added to the 
result to obtain the final answer. 


relative magnitudes are not as stated above, then u = 


The following approximation is used: 
Z 


2k-1 
arctan u 2 Copy U 


C, = .15914533  C, = =,05302625 C, = .03152520 


C, =-.02106178 Cy = .01267292 C,, = --00534860 


Ci3 = .00108423 
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Program Operation 


Program Assembly and Loading 
This program is assembled by the SLAP 2 System (2186651) 
using this symbolic card deck for this program. The SLAP 2 


Loader uses the object card deck, a product of the SLAP 2 
Assembly, to load the program into core memory. 


Program Control 


This program is used as a subroutine by a control program 


Control Program 


Calling Sequence 
Location OP Code Address 


_ Argument (see below) 
+1 TSP ATN (see below) 
+2 Return 
Argument - The fixed point value of the numerator, Y, is in 

the left (L) accumulator and the fixed point 
value of the denominator, X, is in the right (R) 
accumulator. Values Y and X may be at any scale, 
but the scaling must be identical. 


ATN - Symbolic entrance in the arctangent subroutine. 


Output Format 


The value of the arctangent ¥ will be in the left (L) 


accumulator in fractions of a circle scaled at 0. The 
result is not restricted to principal values, but may take 
on any value between -180° and 179.99998°. 





Maximum relative error is less than 107°, 


Approximate running time is 14.5 milliseconds. 
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Memory Requirements 


See Program listing. 
Index registers used - none 


Arctangent - Fixed Point - Card Deck 





The symbolic source card deck for this program is 2186637-1 
Appendices 

Flow Chart 

The flow chart for this program is in Appendix A 

Tables 

Program Listing 

A listing of this program showing memory locations (relative’ 
to zero), octal representation and source card contents (in 


symbolic language format) for each memory location used by 
the program is in Table 1. 
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FLOW CHARTS 
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WAS Y NEGATIVE? 


COMPLEMENT AND STORE ABSOLUTE Y 


WAS X NEGATIVE? 
COMPLEMENT AND STORE ABSOLUTE X i 
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LOAD X, SHIFT RIGHT 1 PLACE 
AND DIVIDE BY Y 


LOAD Y, SHIFT RIGHT 1 PLACE 
AND DIVIDE BY X 


SQUARE ARGUMENT 
AND STORE 


CALCULATE: 7 
arctan X = & Cc. 
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Introduction 


The purpose of the fixed point Log X routine is to calculate 
the logarithm to the base 2, e or 10, of an argument X. 


This mathematical routine has been prepared to fulfill the 
requirements of Specification Number CDV 86638. This drawing, 
together with the reference documents, describes the procedures 
for preparing program card decks and the method of utilizing 
same. 


The program flow chart shown in Appendix A describes the logic 
used to develop the program. 


Mathematical Method 





The argument X is tested and if greater than zero, is 
normalized. This reduces the argument to the form: 


where se XS 1 
and b = binary scaling 
The following approximation is used_to compute log X: 


3 


2K-1 4. 
log, X = Coy U =F 


pe ke NI 
x + Ni/2z 
28853913 
0.96147063 


0.59897865 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS | SIZE | CODE IDENT NO. 


‘ARE IN INCHES AND INCLUDE THICKNESS 638 
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[ARE IN INCHES AND INC 


The result (characteristic) is scaled at 5 and multiplied by 
the proper conversion factor since 


log,A = log, X = log ,2 


The scaling factor is computed and added to the characteristic 
and scaled at 5 for the final answer. 


Program Operation 


Program Assembly and Loading 
This program is assembled by the SLAP 2 system (2186651) using 
the symbolic card deck for this program. The SLAP 2 Loader 


uses the object card deck, a product of the SLAP 2 Assembly, 
to load the program into core memory. 


Program Control 
This program is used as a subroutine by a control program. 
Control Program 


Calling Sequence 


Location OP Code Address 





= Argument (see below) 
2. TSP Entry Points (see below) 
+2 Return 


Argument - The fixed point argument, X, is in the left (L) 
accumulator. The scaling of the argument is in 
the right (R) accumulator scaled at 23. 


Entry Poiats- The program has three entry points. Select one 
of the symbolic entry points listed below 
depending upon the base desired 





Base Entry 


LG2 
LGE 


(CLUDE THICKNESS 


OF PLATING, TOLERANCES ARE: 


Kat 
FRACTIONS = 
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ANGLES = + 1/2° 





Output Format 

The result is in the L accumulator scaled at 5. 

Error Alarm: If the value of the argument is zero or 
negative or if the characteristic is greater than 31,9» 


the overflow indicator will be set and the L and R 
accumulators will contain erroneous data. 


Memory Requirements 


See Program listing. 
Index Registers Used - none 


Log X (2, e, 10) - Fixed Point - Card Deck 





The symbolic source card deck for this program is 2186638-1 


Appendices 


Flow Chart 


The flow chart for this program is in Appendix A. 


Tables 


Program Listing 


A listing of this program showing memory locations (relative 
to zero), octal representation and source card contents (in 
symbolic language format) for each memory location used by 
the program is in Table 1. 
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Introduction 


The purpose of this document is to define the areas common 
to a group of mathematical subroutine written in a floating- 
point format. 


Description 
Floating Point Number 


A floating point binary number (X) is represented by a signed 
proper fraction (M) times some integral power (c) of 2. In 
a normalized floating point number the binary point is posi- 
tioned to the left of the most significant digit of (M). 
Examples: 
X (binary) M (binary) 2° (decimal) 
+001 -100 
-.001 +100 
+100 +100 
1.01 -101 
110. .110 
000. .000 


In the computer a floating point binary number appears in the 
L and R register as follows: 


L-register R-register 





s | Character- | Most significant s | Least significant 
istic part of fraction part of fraction 











23. 22 12 ll 0 23° 22 0 
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The fraction is contained in bits 0 through 11 in the L-register 
and bits 0 through 22 in the R-register. The sign of the 
fraction is contained in the sign bits of both the L and the 
R-registers. The characteristic is formed by adding 2000, 

to the exponent (c). For example an exponent of 100g 
would be represented by a characteristic of 2100g. A zero 
number has no bits in either the L or R-registers. A negative 
number is expressed by complementing the absolute value of the 
fraction and the characteristic together. 


Example: 
L-register R-register 


1 expressed in floating point 20014000 00000000 
-1 expressed in floating point 57764000 40000000 


Subroutine 


Each subroutine is described individually in detail. Associated 
with each subroutine is an identification code. This symbolic 
ID code is used in the calling sequence for entry into the 
subroutine. 


In the calling sequence, the various error returns are 
described. The terms overflow and underflow frequently are 
given as error returns. Overflow refers to the condition in 
which the characteristic exceeded 3777g. Underflow refers to 
the condition in which the characteristic is less than 0. 


The only restrictions for any of the subroutines will be that 
they require that other subroutines be available. For example, 
BCD uses the DVD and FEX subroutines, hence these two sub- 
routines must be available when using the BCD subroutine. 





The floating point subroutines assume all arguments to be in 
the floating point binary format. The output of each of the 
subroutines will be in normalized floating point binary format. 
While the arguments of the subroutines need not be in 
normalized floating point, greater accuracy will result if 
they are normalized. 
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All of the subroutines are written in SLAP-2 (2186651) 
assembly language. 





The subroutines use a restricted set of defined symbols, thus 
allowing them to be assembled together and with other programs. 
The symbols used within any one subroutine are defined by the 
SLAP-2 ID of that subroutine and the symbol "COM." Each 
symbol has as its first three letters the subroutine's 

SLAP-2 ID. For example, the only symbols used in the square 
root subroutine, FSR, are: 





FSR 
FSRA 
FSRB 


FSR9 


The "COM" symbol is used for temporary storage required by 
subroutines common to this system. The requirements of this 
common temporary storage is defined in each subroutine as 
temporary space required in memory. The symbol "COM" appears 
in the location field and the address field of the subroutines 
requiring this temporary storage. 
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To utilize this common storage by a program using more than 
one of these subroutines, the following steps should be taken: 


a) Remove the card, which has the COM symbol defined in the 
location field, from the symbolic source card deck of 
each floating point subroutine. 


Remove the END card from each program. 


Insert one card, which has the COM symbol defined in the 
location field, representing the largest number of 
temporary storage locations required for any one subroutine 


d) Assemble as given under Program Assembly for each 
subroutine. 


System Subroutines 


The subroutines associated with this Floating Point System are 
listed below: 


Program 

Program Number 

Binary Coded Decimal to Floating Point Binary 2186639 
Floating Point Binary to BCD Floating Format 2186640 
Floating Point Add or Subtract 2186641 
Floating Point Multiplication 2186642 
Floating Point Divide 2186643 
Floating Point Square Root 2186644 





Floating Point Sin or Cosine 2186645 
Floating Point Arctangent 2186646 
Floating Point Common Logarithm or Natural Logarithm 2186647 
Floating Point Exponential, E* 2186648 
Floating Point Tangent or Cotangent 2186649 
Double Precision Divide 2186650 
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Identification 


ID: BCD 


Introduction 


The purpose of "Binary Coded Decimal to Floating Point 
Binary" program is to convert a floating point binary 
coded decimal number into binary floating point format. 
This allows floating point numbers to be entered into 
the computer in a decimal floating point notation. 


This program conforms to the specifications of the Floating 
Point System (2186654). 


This mathematical routine has been prepared to fulfill the 
requirements of Specification Number CDV 86639. This drawing, 
together with the reference documents, describes the pro- 
cedures for preparing program card decks and the method of 
utilizing same. 


The program flow chart shown in Appendix A describes the 
logic used to develop the program. 


Mathematical Method 
To change a Binary Coded Decimal floating point number to 


binary floating point number, the following procedure is 
fo Llowed : 


First compute the value of the input mantissa, M,, and the 
input characteristic, C,. 


10 


= 19720 a, 1010t 
Ma i 
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The value of the binary floating point number, X, is 
X= My 10° 


xo My 20a (Log, 10) 


x= M, 2° 24 where Cy is an integer, 


C, 
Using e* subroutine to evaluate 2°4 


200 = M, 2N where Cy is an integer, 


Cc CN 
x=M c 
B M, 2 2 
Combining tems 


X=M 2° whece C=C +¢, 


MMi Me 


N 


Program Operation 


Program Assembly and Loading 


This program is assembled by the SLAP 2 System (2186651) 
using the symbolic card deck for this program and the 
associated subroutines (Section 7.0). The SLAP 2 Loader 
uses the object program card deck, a product of the SLAP 2 
assembly, to load the program into core memory. 





Program Control 


This program is used as a subroutine by a control program. 
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Control Program 
Calling Sequence: 
Location Operation Address 


a TSP BCD (see below) 
atl #5 M,T (see below) 
at2 Return 





BCD - Symbolic entrance in the Binary Coded Decimal to 
Floating Point Binary Subroutine. 


M is the address of the first word of the BCD (Binary Coded 
Decimal) character table of four words which defines one 
floating point number. 


T (04 T4 6) is the index tag modifying M. 


Input Data errors are indicated by an overflow 'ON' condition 
upon returning from subroutine. Overflow will be 'RESET! if 
there are no input data errors. 


The BCD character table is in the following format in memory. 


Location 23 18 17 12 ll 6 5 0 
M 














Sign of the mantissa in BCD. It can be +, -, 
or space. 


Most significant digit of the mantissa in BCD 


Least significant digit of the mantissa in BCD. 


Sign of the characteristic in BCD. It can be +, 
-, or space. 


Most significant digit of the characteristic in BCD. 
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Output Format 


Normal return: 


The normalized floating point number equivalent to the speci- 
fied BCD number will be in the L and R registers upon return. 


Error Return: 
a. The characteristic of the BCD number is greater than or 
equal to 1024. The L register will be positive upon 


return. 


The characteristic of the BCD number is less than -1024. 
The L register will be negative upon return. 


Memory Requirements 
Program: See program listing 

Temporary Space: 14 core memory locations 

Restrictions 

This program uses the following floating point subroutine: 


Name 1D Number 


Exponential, E* 2186648 


Timing 





Approximate timing: 100.9 milliseconds 


“Binary Coded Decimal to Floating Point Binary" Card Deck 





The symbolic source card deck for this program is 2186639-1 . 
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Flow Chart 


Flow Chart for this program is in Appendix A. 







Tables 


Program Listing 






A listing of this program showing memory locations (relative 
to zero), octal representation and source card content (in 
symbolic language format) for each memory location used by 
the program is in Table l. 
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Identification 


ID: BIN 


Introduction 


The purpose of the Floating Point Binary to BCD Floating 
Format program is to convert a floating point binary quantity X 
into a Binary Coded Decimal floating point format. 


This program conforms to the specifications of the Floating 
Point System (2186654). 


This mathematical routine has been prepared to fulfill the 
requirements of Specification Number CDV 86640. This drawing, 
together with the reference documents, describes the procedures 
for preparing program card decks and the method of utilizing 
same. 


The program flow charts shown in Appendix A describes the 
logic used to develop the program. 


Mathematical Method 


To convert a binary floating point number to a Binary Coded 
Decimal floating point format, the following procedure is 
followed. 


The floating point binary number is represented thusly: 


i =~ m2? 
x = mio? 108 2 
x = M1022 1073 where Dj is an integer 


D3 is a fraction 
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Using the EX subroutine to evaluate 1023 = e231n10 






X = mig22_P3In10 


x= F 1072 where F = MeD31"10 


F< 10 







1s 





IfF = 1, set F= F + 10 and D) = Dy 






Since ib 


F=)) x10! 


i=l 










successive multiplication by 10 yields the M, terms. 







Since 


Bei 
Do = Do 6,10 
ii 









successive division by 10 yields the C, terms. 


Program Operation 
Program Assembly and Loading 


This program is assembled by the SLAP 2 System (2186651) 
using the symbolic card deck for this program and the 
required associated subroutine (Section 7.0). The 

SLAP 2 Loader uses the object program card deck, a product 
of the SLAP 2 Assembly, to load the program into core 
memory. 











Program Control 


This program is used as a subroutine by a control program. 
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Control Program 


Calling Sequence: 


With the floating point binary quantity, X, in the left (L) 
and right (R) accumulators 





Location Operation Address 




















TSP BIN (see below) 
+1 KS M,T (see below) 
+2 Return 





BIN - Symbolic entrance in the Floating Point Binary to BCD 
Floating Format subroutine. 


M is the address of the first word of the BCD (Binary Coded 
Decimal) character table of four words which is defined by 
the given floating point binary number. 











T (O£ T€ 6) is the index tag modifying M. 


Overflow 'ON' indicates input error. Overflow will be reset 
on the normal return. 





The BCD character table is in the following format in memory: 


















23 





location 18 O67 2 








M 
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Sign of the mantissa in BCD. It can be - or space. 


Most significant digit of the mantissa in BCD. 


- Least significant digit of the mantissa in BCD. 


Sign of the characteristic in BCD, It can be - or 
space 
Most significant digit of the characteristic in BCD 


Least significant digit of the characteristic in BCD 


Output Format 

The binary Coded Decimal floating point format of the given 
binary floating point number has been stored in table 
specified by the address M of the calling sequence. 


The contents of the (L) and (R) registers have been destroyed 
by the subroutine. 


Memory Requirements 
Program: See Program Listing 
Temporary Space: 12 core memory locations 


Restrictions 





This program uses the following floating point subroutine: 


Name 1D Number 


Exponential, EX 2186648 
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Timing 
Approximate Timing: 99.6 milliseconds 


Floating Point Binary to BCD Floating Format Card Deck 





The symbolic source card deck for this program is 2186640-1 


Appendices 


Flow Chart 


Flow Chart for this program is in Appendix A. 


Tables 


Program Listing 


A listing of this program showing memory locations (relative 
to zero), octal representation and source card contents (in 
symbolic language format) for each memory location used by the 
program is in Table 1. 
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1S THE NUMBER = 0 7 





18 THE NUMBER 
NEGATIVE 7 








=D, + Dy 
COMPLEMENT Dy IS INT., 
THE NUMBER a Dg 38 FRACT. 


ea Fae 


SET THE SIGN D4 ® Dg In 10 
OF THE SCALED AT 
MANTISSA, Two 


Mg, MINUS 
[NORMALIZE D, 


GET THE ‘ LET ya NUMBE! 
MANTISSA, M lor SHIFTS TO 


a” [NORMALIZE 


FLOAT NORMALIZED 
GET THE COMPLEMENT Dy. (102G-y)= CHAR. 
CHARACTERISTIC ID Dy = MANT. 


q 1 


ueustHe | SET SIGN OF 
CHARACTERISTIC CHAR. Cs, | 
|D = CHAR. - 1024 NEGATIVE 







































































D LESS THAN ZERO ? aes 


HSE "|a|49671| 2/86640 


rucrons == _mousse2ye [SAE [SET 7 
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WAS Cy = NEGATIVE 


y= INTEGER 
PART OF (10f) 
FRACTIONAL 





STORE SIGN OF 
CHARACTERISTI 


‘SIGN 7 
YES 








COMPLEMENT, 











STORE XIN 
Derr. MOST SIGNIF. 


COUNTER BY PART OF 
ONE. CHAR. 


er! 


18 DIGIT COUNT= 107 


INCREASE 

















2 
(suieT THE 
CHAR. DIGITS 
RIGHT 6 
PLACES 


STORE SIGN OF 


IMANTISSA WITH 
JMANTISSA DIGITS 


1S CHARACTERISTIC 
ZERO? 


SET SIGN OF 
CHAR, PLUS 





INCREASE 
DIGIT 
COUNTER BY 
ONE 























DIGIT COUNT = 3 ? 








STORE FLOATE! 
NUMBER IN 
SPECIFIED 
ADDRESS 





RADIO CORPORATION OF AMERICA 


Frameamesne”” |Al49671| -/B6OCFO 
a a 7 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS ‘a | ‘CODE IDENT NO. 


Met 
FRACTIONS = 
DEL 12633A (2/65) 


OK = ot 
ANGLES = 172° 























TABLE 1 


SATURN GROUND COMPUTER SYSTEM 


FLOATING POINT BINARY TO 


BCD FLOATING FORMAT 


PROGRAM LISTING 


RADIO CORPORATION OF AMERICA 









CODE IDENT NO. 4 


49671| 2/86640 | 
a a Ce 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS. 
ARE IN INCHES AND INCLUDE THICKNESS 
(OF PLATING, TOLERANCES ARE: 


Kat 2K 
FRACTIONS = ANGLES = 1/22 
DEL 1263-30 (2/65) 

























(Sw) (52) ‘T2- ay GaTv9s OF NI a) 


($0) Gio *T2- ay gatas OF NT (80) 


iS) 
+ 
‘0 
‘0 
XQ 
N 


lava SH "(2 90700 


‘or uy aaqvas (2 90790 
‘yaya ANVOLSINOTS 45¥37 “(2 90700 

*IN3N@GX3 3801S 
S1¥319¥avHD 13s 
*INBNBGX3 3AL1v93N LN3W31GHOO 





CNBIS OTLS1¥aLOVAWKO BHU 13s 
"OU1S1H3L9vENHO 3ML S¥TENN 


CODE IDENT NO. 
49671 


"OuY 3HL MOMS YSSTLNVW 3H, 3A0NaY 


'0 Ly WSSILNYM ams 37¥9S 
*OLISIaL9¥eNHO BML ind LovaLea 


x3 3HL OL Ba4SNVL 


“Sauon 028 KI 3002 8uaz a0v14 
*ouaz SI (N3WndwY 3KL 
sauaz St ANaKNOBY BML 41 93ND 


° . 
098 01 AYYNIG INIod ONTLYOT 





woe 
adat| wwi00 | -907 
ANBLNOD Q¥YD 91708NAS sau | awonam | isn 





UNLESS OTHERWISE SPECIFIED: DIMENSIONS | 


[ARE IN INCHES AND INCLUDE THICKNESS 
(OF PLATING. TOLERANCES ARE: 


‘VORB 30 NOLLVHOSHOO O10Ve | 


= 1/22 


ANGLES: 





DEL 1263-3A (2/65) 


ty Av Wed ‘1ondeud 7103 
(sio (si) ty LY WY Londada WI iave 


“yoLv wy L9ndeud Wiiavd 


(S60 ($2) *6T= LY WY aungnoD 


(S81) (SW) *6T> LY WY aungwa2 


‘avy 39vos A73LNI8sav 

“0 ayy 33¥98 
*tt ay aaqwas sI_y 
*V 48 YSSIINYH 3HL 139 





“YSSLLNYH 
any ¥ 30 





TNT SGedsv 38g1S 
“BNIunewans TILNaNBaxa 


INawa7gH0o 


ANY 3HL 3003m 


41WS_9NIZ1- THN) -920T 
VAYHD 3ML aLndWad 
1vO1s 81 3¥vaaua 
OF a (sa) atv9s, 


"Or NT (20) aZITWHUON 


(sw) (Sw) '2 LY ggqwas Ot WI (20) 
*kongeus Wliuvd °2 1¥ gaqvds Ot NT (ed) 


zoe0 0 cor 


stexnle 


ate 
gexnle 


'YORIaMIY 30 NOLLVOSUC O1GVE 


al aser!| 2/8 6640 


‘SHEET /5 


ANGLES = #1/2° 


=s 





FRACTIONS 





“SANIW ST_NOIS 41 3aMMi GyBM vied OINI 
‘D1istwaiovaWKd 3M, 40 NDIS 30URH 
“ouaz SI_iNaNedxa 

1 S016 SI IN3NGGKA 48 MOIS 


*20VaS YATTON SI ROIS 
“SAWIM SINOIS 41 3N@ QUBM Yivd @INT 

NVWaML 48 NDIS 208aM 
1d xan was NUNLZE *ON 
saassapeue naae “3am S11910. 43) 41 393K 





“ANG AB LN 41010 3KL 
“e OMBK 
“2 Oyen 
“Saovie 9 1437 33uML Ons SdUBK 





*GMon GuIML 3Hi BLNI L101 aHL 30uaH 
*ouaz st i1910 41 x9aH9 


Quen viva GMIHL SNL 48 waLovHWHD 
ANVOLSINOIS 1S¥37 BML NI LTOTO aL 39vIG 


souen vivo a98 
BML NI 4437 9 SuaLovuNHD LaTHS 

‘yay sor 

"407 48 STYH ANVOTSINDIS 1S¥37 31ndHaD 


+401 48 40H ANVOIGINOIS aH, 3uBLS 
‘wiv gatos 4 oF 


*O Av aav9s 38 O41 4 LatKS 
“wy LV NOLLOVES Ine ¥SYH 
*AwaNoax3 TWI930 3M1 01 aNB GOV 


“BNB_NYHL S$37 SI 4 31 yasswyaL 





OTeXNIB 


'YORISMV 30 NOLLVEOSUOD O1GVE 


‘a as671| 2/86640 


SCALE 


ANGLES = + 1/2° 





DEL 1263-34 (2/65) 


oT NI (ea) MBs wordy, 
OF NT (80) wed Yardy OF 


lave S17 
“Sule 2% yanay us 


iuva $7 


‘ayen 098 ayaz 


*SdWOH ViVO Muungs ONY GaIM, aHL 3uBLS 


*SdWOM Viva GNBOIS GNY IswI4 aHL 3uBIS 


Suen viva » aML auBLS 
*aBWHL SI 1Nn@D 41910 41 ¥93HO 


*3NB AB ANned 11910 3K ane 


*A3NBdX3__BKL 
“duon viva Hiunes 3HL 3z1 
*eu3z inne2 110) 








gorceezzo 
zszs90680 


0 20 
‘0ozozozo 
‘ozozozczo 
ez0z0z000 


weexne 


fore) 
senda 
rer) 
2-409 
s2exnle 


20020000 
0000000 
90909000 
soccerze 
zseaooes 
“ 0 


ry T 
0000002 
oo00ccce 
44220000 





9920 0 10y 


‘WOREaRNY 30 NOLLVUOSUCO O1OVE 


Z/8664O0 | 
ser ya 


‘CODE IDENT NO. 


A| 49671 





DEL 1263-3A (2/65) 


x34 


00000 0 00 


'VORIINY 40 NOLLVEOSUOO O1VE 


2/86640 | 


‘CODE IDENT NO. 
49671 


‘SIZE 


5 
G 
5 
4 
§ 


21/20 





ANGLES: 


=+ 


FRACTIONS 
DEL 1263-34 (2765) 


SECTION 9 
FLOATING POINT ADD OR (SUBTRACT) 


REVISIONS 
DESCRIPTION DATE | APPROVED 











FIRST MADE FOR 


NEXT ASSY USED ON 














CONTRACT; NO RADIO CORPORATION OF AMERICA 
NAS 8-13007 NEW YORK, NY—VAN NUYS PLANT 
DRAWN DATE 
ne 2.8€p.65| SATURN GROUND COMPUTER SYSTEM 
EME ATE MATHEMATICAL ROUTINE 
AC whe... B2ceus| FLOATING POINT ADD OR (SUBTRACT) 


DESIGN ACTIVITY APPD . | SIZE | CODE IDENT NO. 
6 2/8 664/ 


JH Vandufed F Saerbs| 
SHEET 1 OF 








DEL 1263-1 (2-65) 


LIST OF MATERIALS OR PARTS LIST 
QTY REQD PART OR * VENDOR ITEM — SEE SOURCE CONTROL OR SPEC CONTROL DWG 
IDENTIFYING NO. NOMENCLATURE OR DESCRIPTION ‘SPECIFICATION 


2186653 SATURN GROUND COMPUTER SYSTEM 
[SPECIFICATION FOR SATURN V 
COMPUTER PROGRAM 

2186654 SATURN GROUND COMPUTER SYSTEM 
MATHEMATICAL ROUTINE 

FLOATING POINT SYSTEM 

2186651 SATURN GROUND COMPUTER SYSTEM 
SLAP-2 SYSTEM 

2186641-1 FLOATING POINT ADD OR (SUBTRACT) 
SYMBOLIC CARD DECK 

CDV 86641 CEI DETAIL SPECIFICATION 
MATHEMATICAL ROUTINE. 

FLOATING POINT ADD OR (SUBTRACT) 
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Identification 


ID: FAD or (FSB) 
‘ Introduction 


The purpose of the Floating Point Add or (Subtract) program 
is to compute the floating point sum, X + Y, or the 
remainder, X ~ Y, where both X and Y are floating point 
arguments. 


This program conforms to the specifications of the Floating 
Point System, 2186654. 


This mathematical routine has been prepared to fulfill the 
requirements of Specification Number CDV 86641. This drawing, 
together with the reference documents, describes the pro- 
cedures for preparing program card decks and the method 

of utilizing same. 


The program flow chart shown in Appendix A describes the 
logic used to develop the program. 





Program Operation 


Program Assembly and Loading 


This program is assembled by the SLAP 2 System (2186651) 
using the symbolic card deck for this program. The SLAP 2 
Loader uses the object program card deck, to load the 
program into core memory. 





Program Control 


This program is used as a subroutine by a control program. 
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Control Program 


Calling Sequenc: 





With the floating point argument, X, in the left (L) and 
right (R) accumulators, 


Location Operation Address 


a TSP FAD or (FSB) 
(see below) 
5 M, T (see below) 


avy 
at2 Return 
FAD - Synbolic entrance in the Floating Point Add subroutine 


FSB - Symbolic entrance in the Floating Point Subtract 
subroutine. 


M is the address of Y the addend or (subtrahend) 
T (0$ T< 6) is the index tag modifying M 
Output Format 

Error return: 

An error return is indicated by overflow 'ON'. 


a. Overflow condition exists. The (L) register will be 
positive upon return, 


Under flow condition exists. The (L) register will be 
negative upon return, 


Normal return: 
Normal condition is indicated by overflow 'OFF'. 
The sum or (remainder) of the two arguments, X and Y, will 


be in the L and R registers upon return, The result will 
be in normalized floating point format. 


ARE IN INCHES AND INCLUDE THICKNESS 
TOLERANCES ARE: 











Memory Requirements 


Program: see program listing 


Temporary Space: 7 core memory locations 


Restrictions 


None 
Timing 
Approximate Timing - 5.25 milliseconds 


Floating Point Add or (Subtract)Card Deck 





Symbolic source card deck for this program is 2186641-1 
Appendices 

Flow Chart 

Flow chart for this program is in Appendix A. 

Tables 

Program Listing 

A listing of this program showing memory locations (relative 
to zero), octal representation and source card contents 


(in symbolic language format) for each memory location 
used by the program is in Table 1. 
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COMMAND FOR CMB IF FSB 
Y ARGUMENT 








(oeT THEY 
18 X NEGATIVE? IMANTISSA, M 
JAND SCALE 
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ae 


N= CX- CY 


PICKUP Y 
ARGUMENT 





GET THE X 

CHAR Cx 18 Y NEGATIVE YES 

AND SCALE 

AT 23 YES wnze? 
NO 








|GET X MANTISSA,| jCOMPLEMENT 
MX AND SCALE Y 


Rr ONE COMPLEMENT 














BUILD LDB 
|COMMAND TO 
PICK UP ¥ 
ARGUMENT 


GET THEY 
CIIAR., CY. 
SCALE CY 
AT 23 
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IS Z EQUAL TO ZERO ? 


18 N NEGATIVE? 








a2 

BUILD SBA 
COMMAND TO 
SHIFT N 
PLACES 
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13 N>457 
































7 DERFLOW. ERROR 
18 CHAR, NEGATIVE? BET L REGISTER) RETURN 
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BCuAR 2 s0ie) er a necisren 


COMBINE Z 
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CHAR 
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OEE NO, RADIO CORPORATION OF AMERICA 
NAS 8-13007 NEW YORK, NY—VAN NUYS PLANT 
DRAWN DATE 
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CHECKED DATE 
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DESIGN ACTIVITY APPD el | SIZE | CODE IDENT NO. 
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LIST _OF MATERIALS OR PARTS LIST. 
QTY REQD’ PART OR * VENDOR ITEM — SEE SOURCE CONTROL OR SPEC CONTROL DWG. 
IDENTIFYING NO. NOMENCLATURE OR DESCRIPTION SPECIFICATION 


2186653 SATURN GROUND COMPUTER SYSTEM 
SPECIFICATION FOR SATURN V 
COMPUTER PROGRAM 

2186654 SATURN GROUND COMPUTER SYSTEM 


MATHEMATICAL ROUTINE 
FLOATING POINT SYSTEM 


2186651 SATURN GROUND COMPUTER SYSTEM 
SLAP-2 SYSTEM 

2186642-1 FLOATING POINT MULTIPLICATION 
SYMBOLIC CARD DECK 

CDV 86642 CEI DETAIL SPECIFICATION 
MATHEMATICAL ROUTINE 
FLOATING POINT MULTIPLICATION 
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Identification 


ID: FMP 


Introduction 


The purpose of the Floating Point Multiplication program is 
to compute the floating point product XY where X and Y are 
both floating point numbers. 


This program conforms to the specifications of the Floating 
Point System, 2186654. 


This mathematical routine has been prepared to fulfill the 
requirements of Specification Number CDV 86642. This 
drawing, together with the reference documents, describes 
the procedures for preparing program card decks and the 
method of utilizing same. 


The program flow chart shown in Appendix A describes the 
logic used to develop the program. 


Mathematical Method 
For ease of notation, let M designate the most significant 
part and L the least significant part of the multiplicand X. 
Similarly, let R designate the most significant part and 
S$ the least significant part of the multiplier Y. 

XY = MR + (MS + LR) 
Program Operation 


Program Assembly and Loader 


This program is assembled by the SLAP 2 System (2186651) 
using the symbolic card deck for this program. The SLAP 2 
Loader uses the object program card deck, to load the 
program into core memory. 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS | SIZE | CODE IDENT NO. 
ARE IN INCHES AND INCLUDE THICKNESS 
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Program Control 


This program is used as a subroutine by a control program. 


Control Program 
Calling Sequence: 
With the multiplicand X in the left (L) and right (R) 
accumulators, 

Location Operation Address 


a TSP FMP (see below) 
if 5 M,T (see below) 
2 Return 


art 
a+ 


FMP - Symbolic entrance in the Floating Point Multiplication 
subroutine. 


M is the address of the multiplier, Y, 
T (0S TS 6) is the index tag modifying M. 


Output Format 


Error return: 
An error return is indicated by overflow 'ON'. 


a.  Underflow condition exists. The (L) register will be 
negative upon return. 


b. Overflow condition exists. The (L) register will be 
positive upon return. 

Normal return: 

Normal return is indicated by overflow 'OFF'. 


The product of XY will be in the (L) and (R) registers 
upon return. The product will be in normalized floating 
point format. 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS | SIZE | CODE IDENT NO. 
{ARE IN INCHES AND INCLUDE THICKNESS 2/86 Par h ae 
OF PLATING. TOLERANCES ARE: 








Memory Requirements 
Program: see program listing 
Temporary Space: 10 core memory locations 
Restrictions 
None 
fet 


Approximate timing: 6.4 milliseconds 





Floating Point Multiplication Card Deck 





Symbolic source card deck for this program is 2186642-1. 
Appendices 

Flow Chart 

Flow chart for this program is in Appendix A. 

Tables 

Program Listing 

A listing of this program showing memory locations (relative 
to zero), octal representation and source card contents (in 


symbolic language format) for each memory location used by 
the program is in Table 1. 
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LIST OF MATERIALS OR PARTS LIST 
QTY REQD PART OR % VENDOR ITEM — SEE SOURCE CONTROL OR SPEC CONTROL DWG 
[503 [502) IDENTIFYING NO. NOMENCLATURE OR DESCRIPTION ‘SPECIFICATION 
2186653 SATURN GROUND COMPUTER SYSTEM | 
SPECIFICATION FOR SATURN V 
COMPUTER PROGRAM 
2186654 “SATURN GROUND COMPUTER sysTEM | 
MATHEMATICAL ROUTINE 
FLOATING POINT SYSTEM 
2186651 SATURN GROUND COMPUTER SYSTEM 
SLAP-2 SYSTEM 
2186643-1 FLOATING POINT DIVIDE — 
SYMBOLIC CARD DECK 
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MATHEMATICAL ROUTINE 


FLOATING POINT DIVIDE 


Prepared Under Contract NAS 8-13007 
for 
National Aeronautics and Space Administration 
Marshall Space Flight Center 
Huntsville, Alabama 
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Identification 


ID: FDV 
Introduction 


The purpose of the Floating Point Divide program is to 
compute the floating point quotient of X + Y where both 
X and ¥ are floating point numbers. 


The program conforms to the specifications of the Floating 
Point System, 2186654. 


This mathematical routine has been prepared to fulfill the 
requirements of Specification Number CDV 86643. This 
drawing, together with the reference documents, describes 
the procedures for preparing program card decks and the 
method of utilizing same. 


The program flow chart shown in Appendix A describes the 
logic used to develop the program. 


Mathematical Method 


For ease of notation, let M designate the most significant 
part and L the least significant part of the divisor Y. 


tb AO] 


Program Operation 





Program Assembly and Loader 


This program is assembled by the SLAP 2 System (2186651) 
using the symbolic card deck for this program. The SLAP 2 
Loader uses the object program card deck to load the 
program into core memory. 
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UNLESS OTHERWISE SPECIFIED: DIMENSIONS | SIZE | CODE IDENT NO. 
[ARE IN INCHES AND INCLUDE THICKNESS 
‘OF PLATING, TOLERANCES ARE: 


Kt 20 = 
FRACTIONS == ANGLES=++1/2° 
DEL 1263-3A (2/65) 





Program Control 


This program is usec as a subroutine by a control program. 
Control Program 
Calling Sequence: 


With the dividend X in the left (L) and right (R) accumulators 


Location Operation Address 


a TSP FDV (see below) 
atl eS, M,T (see below) 
at+2 


Return 


FDV - Symbolic entrance in the Floating Point Divide 
subroutine. 


M is the address of the floating point divisor Y. 
T (O€ T¢ 6) is the index tag modifying M. 


Output Format 


Error return: 


An error return is indicated by over*low 'ON'. 


The divisor is zero. The L and R registers are both 
zero upon return. 


Underflow condition exists. The L register will be 
negative upon return. 








Overflow condition exists. The L register will be 
positive non-zero upon return. 


Normal return: 
A normal return is indicated by overflow 'OFF'. 
The quotient of X + ¥ will be in the L and R registers upon 


return. The quotient will be in normalized floating point 
format. 
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Memory Requirements 

Program: see program listing 

Temporary Space: 10 core memory locations 
Restrictions 

None 

Timing 


Approximate timing: 9.85 milliseconds 





Floating Point Divide Card Deck 

Symbolic source card deck for this program is 2186643-1 
Appendices 

Flow Chart 

Flow chart for this program is in Appendix A. 

Tables 

Program Listing 

A listing of this program showing memory locations (relative 
to zero), octal representation and source card contents (in 


symbolic language format) for each memory location used by 
the program is in Table 1. 
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504] 


503/502 a IDENTIFYING NO. 
2186653 


%* VENDOR ITEM — SEE SOURCE CONTROL OR SPEC CONTROL DWG 





NOMENCLATURE OR DESCRIPTION 
SATURN GROUND COMPUTER SYSTEM 


SPECIFICATION 





SPECIFICATION FOR SATURN V 





COMPUTER PROGRAM 





2186654 


SATURN GROUND COMPUTER SYSTEM 





MATHEMATICAL ROUTINE 





FLOATING POINT SYSTEM 





2186651 


SATURN GROUND COMPUTER SYSTEM 





SLAP-2 SYSTEM 





2186644-1 
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Identification 


ID: FSR 


Introduction 


The purpose of the Floating Point Square Root program is to 
compute the floating point square root of X where X is in the 
floating point format. 


The program conforms to the specifications of tie Floating 
Point System 2186654. 


This mathematical routine has been prepared to fulfill the 
requirements of Specification Number CDV 86644. This drawing, 
together with the reference documents, describes te procedures 
for preparing program card decks and the method of utilizing 
same. 


The program flow chart shown in Appendix A descri‘es the logic 
used to develop the program. 


Mathematical Method 


Assume that the argument X is normalized although this is not 
a program requirement. If the argument is not normalized, it 
will immediately be normalized. 


Unpack X, 


xh/2 2 gyt/? ell? where 1> € > 125 





=a 


The firet fractional approximation, yqs to £!/? is given by 
£+ 1/4 - 1/2(£ - 1/2)%. This is then followed by three 
Newtonian iterations of the form: 
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Yep tle, oy 
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Program Operation 

Program Assembly and Loader 

This program is assembled by the SLAP 2 System (2186651) 
using the symbolic card deck for this program. The SLAP 2 


Loader uses the object program card deck, to load the program 
into core memory. 


Program Control 


This program is used as a subroutine by a control program. 


Control Program 


Calling Sequence: 


With the argument X in the left (L) and right (R) 
accumulators, 


Location Operation Address 





a TSP FSR (see below) 
atl Return 


FSR - Symbolic entrance in the Floating Point Square 
Root subroutine. 


5.0 Output Format 


Error Return: 
An error return is indicated by overflow 'ON'. 


The argument X was negative. The L and R register will 
contain the negative argument. 


Normal Return: 
A normal return is indicated by overflow 'OFF'. 


The square root of the argument X will be in the L and R 
registers upon return. The square root will be in normalized 
floating point format. 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS | SIZE | CODE IDENT NO. 
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Memory Requirements 


Program: see program listing. 
Temporary Space: 12 core memory locations. 
Restrictions 

None 

Timing 

Approximate timing: 26.6 milliseconds. 


Floating Point Square Root Card Deck 





Symbolic source card deck for this program is 2186644-1. 
Appendices 

Flow Chart 

Flow chart for this program is in Appendix A. 

Tables 

Program Listing 

A listing of this program shown memory locations (relative 
to zero), octal representation and source card contents (in 


symbolic language format) for each memory location used by 
the program is in Table 1. 
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SATURN GROUND COMPUTER SYSTEM 
FLOATING POINT SQUARE ROOT 
PROGRAM LISTING 
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Identification 


ID: FSN or (FCS) 


Introduction 


The purpose of the Floating Point Sine or (Cosine) is to 
compute the floating point sine or (cosine) of a floating 
point radian argument X. 


The program conforms to the specifications of the Floating 
Point System 2186654. 

This mathematical routine has been prepared to fulfill the 
requirements of Specification Number CDV 86645. This 
drawing, together with the reference documents, describes 
the procedures for preparing program card decks and the 
method of utilizing same. 


The program flow chart shown in Appendix A describes the 
logic used to develop the program. 


Mathematical Method 
To evaluate sin X, the argument X is first adjusted in the 


following manner: For positive X, change from radians to 
parts of a circle. 


& = K + L where K is an integer and L < 1. 
4 L= 1+ F where I is an integer and F< 1. 


If1l=0, let z= 
If1= 1, let Zz 
If1l= 2, let Zz 


If1 = 3, let Z 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS | SIZE | CODE IDENT NO. | 
‘ARE IN INCHES AND INCLUDE THICKNESS 


OF PLATING. TOLERANCES ARE: A| 49671 


OK ot 


ANGLES=+1/2° [SCALE 













Mathematical Method (Cont'd) 


Sin X = ggat 


nei 



























0.999999999979082 
-0. 166666666092171 
00083333307 30723 
0.00019 8408338222 
0..000002752401177 
-0.000000023868930 
To compute cosin X, the following relationship is used: 
cosin X = sin (x + 5] 
For negative X, the following relationships are used: 
sin (+ X) = - sin X 
Sin (- X) = cos X 
Program Operation 


Program Assembly and Loader 


This program is assembled by the SLAP 2 System (2186651) 
using the symbolic card deck for this program. The SLAP 2 
Loader uses the object program card deck to load the 

program into core memory. 





Program Control 


This program is used as a subroutine by a control program. 
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Control Program 


Calling Sequence: 
With the angle in the left (L) and right (R) accumulators, 
Location Operation Address 


a TSP FSN or FCS (see below) 
atl Return 


FSN - Symbolic entrance in the Floating Point Sin subroutine. 
FCS - Symbolic entrance in the Floating Point Cosin subroutine. 
Output Format 

Error Return: 

An error return is indicated by overflow 'ON'. 

the |x| >2™*. 

Normal Return: 

A normal return is indicated by overflow 'OFF'. 

The sine or (cosine) of the argument X will be in the L and 

R registers upon return. The sine or (cosine) is in 


normalized floating point format. 


Memory Requirements 


Program: see program listing 





Temporary space: 12 core memory locations 
Restrictions 

None 

Timing 


Approximate timing: 39.3 milliseconds 
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Floating Point Sine or (Cosine) Card Deck 





Symbolic source card deck for this program is 2186645-1. 


Appendices 


Flow Chart 


Flow chart for this program is in Appendix A. 


Tables 


Program Listing 


A listing of this program showing memory locations (relative 
to zero), octal representation and source card contents (in 
symbolic language format) for each memory location used by 
the program is in Table 1. 
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AT TWO abi) 


SCALE MOST 
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OF MAT TWO 


REMOVE BITs 21 
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BITS). CALL THIS 
FRACTION F 
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RADIANS, ADJUST f COMPUTE THE 
R=f(3) SCALED cr CHARACTERISTIC 
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PICK UP COMPUTE MERGE THE 
QUADRANT ef CHAR. AND 
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YES NO 





18 THIS SIN? 













ADD 1 (90 DEGREES) 
TO QUADRANT 
BITS 








MAKE SIN X 
NEGATIVE. 

















QUADRANT 3 OR 4? 





COMPUTE 
& pint 
oxic es 8! 







SET SIGN 


INDICATION in p 
NEGATIVE [SCALED ATONE | 
























NORMALIZE 

SIN X, 
SET SIGN LET S BE NO. 
INDICATION SHIFTS REQUIRED| 









POSITIVE 

























QUADRANT 2 AND 4? 
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FLOATING POINT SYSTEM 
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SLAP-2 SYSTEM. 
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FLOATING POINT ARCTANGENT 





SYMBOLIC CARD DECK 
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Identification 

ID: FRT 

Introduction 

The purpose of the Floating Point Arctangent is to compute 
the floating point arctangent of the quotient (X + ¥). X and 


Y¥ are both floating point arguments of the subroutine. 


The program conforms to the specifications of the Floating 
Point System 2186654. 


This mathematical routine has been prepared to fulfill the 


requirements of Specification Number CDV 86646. This drawing, 
together with the reference documents, describes the procedures 
for preparing program card decks and the method of utilizing 
same. 


The program flow chart shown in Appendix A describes the logic 
used to develop the program. 


Mathematical Method 


x= 2%. x where 1> pl 





y= 2". y where 1 > |y| 


forN>M, Z= 5 and K= 0. 


NeM, Z= and K= 1, 





for | <(v2-1), t=ZandL= 





kl Scvr-n, re GR and L= 1, 
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Mathematical Method (Cont'd) 


Compute 1, 





do+t? (a1 + (2 2K) +x ( 





2 2 
th+d, tet (to +d 


3), 
0.999999999996107 
-0.015585371018178 

2100554087 165198 
1.621023833634443 

0.585315 135071831 
-0.419003002282544 


Final adjustment is made for which quadrant was specified by 
X and Y as follows: 


for Y>0, X>0, arctan 
Y>0, X<0, arctan 
¥<0, X<0, arctan 


Y<0, X>0, arctan 


ee Ted ed 


Program Operation 
Program Assembly and Loader 


This program is assembled by the SLAP 2 System (2186651) 
using the symbolic card deck for this program and the 
associated subroutine (Section 7.0). The SLAP 2 Loader 
uses the object program card deck, to load the program 
into core memory. 
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Program Control 


This program is used as a subroutine by a control program. 


Control Program 


Calling Sequence: 
With the floating point argument ¥ in the left (L) and 
right (R) accumulators: 

Location Operation Address 


TSP FRT (see below) 
#5, M,T (see below) 
Return 





+1 
+2 


FRT - symbolic entrance in the Floating Point Arctangent 
subroutine. 


M is the address of the argument X 

T (02 Té 6) is the index tag modifying M. 
Qutput Format 

Error Return: 

An error return is indicated by overflow 'ON'. 


Both arguments equal zero X= 0, Y= 0. The L register is 
Zero upon return. 


Normal Return: 





A normal return is indicated by overflow 'OFF'. 


The arctangent of ¥ = X will be in L and R registers upon 
return. The arctangent will be in normalized floating point 
radians. The range of the arctangent is 0 to 2m radians. 
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Memory Requirements 


Program: see program listing 


Temporary Space: 14 core memory locations 


Restrictions 


This program uses the following floating point subroutine: 


Name D Number 
Double Precision Divide DvD 2186650 


Timing 
Approximate timing: 49.2 milliseconds maximum. 


Floating Point Arctangent Card Deck 





Symbolic source card deck for this program is 2186646-1 
Appendices 

Flow Chart 

Flow chart for this program is in Appendix A. 

Tables 

Program Listing 

A listing of this program showing memory locations (relative 
to zero), octal representation and source card contents 


(in symbolic language format) for each memory location used 
by the program is in Table 1. 
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Identification 


ID: FLG or (FLN) 


Introduction 


The purpose of the Floating Point Common Logarithm or 
(Natural Logarithm) is to compute the floating point common 
log or (natural log) of a number X which is in floating 
point format. 


The program conforms to the specifications of the Floating 
Point System 2186654. 


This mathematical routine has been prepared to fulfill the 
requirements of Specification Number CDV 86647. This drawing, 
together with the reference documents, describes the pro- 
cedures for preparing program card decks and the method of 
utilizing same. 


The program flow chart shown in Appendix A describes the 
logic used to develop the program. 


Mathematical Method 


The argument X is reduced as follows to yield the natural 
log of X. 


x= 2. F where 1/22 FC 1 
In x = N 1n2 + InF 


InF= in Vi/2+1n (3) where T = 





i-f 


w(tet)er fase? (0+ ep 
(+ T°) 


A = 1.999999999491255 
B = 0.109078890502997 
C = -0.777314001005492 
D = -1.394065145176107 
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Mathematical Method (Cont'd) 


To compute the common log, the following relationship is used: 
log X = log e InX 

Program Operation 

Program Assembly and Loader 


This program is assembled by the SLAP 2 System (2186651) 
using the symbolic card deck for this program and the 
associated subroutine (Section 7.0). The SLAP 2 Loader 
uses the object program card deck, a product of the SLAP 2 
Assembly, to load the program into core memory. 


Program Control 


This program is used as a subroutine by a control program. 


Control Program 


Calling Sequence: 


With the floating point argument X in the left (L) and right 
(R) accumulators, 


Location Operation Address 


a TSP FLG or (FLN) 
(see below) 
atl Return 





FLG ~ symbolic entrance in the Floating Point common 
logarithm subroutine 


FIN - symbolic entrance in the Floating Point natural 
logarithm subroutine 
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Output Format 


Error Return: 


An error return is indicated by overflow 'ON'. 


a. The argument X is zero. The L register contains 
(40000000) and the R register contains (40000001) 
upon return. 

The argument X is negative. The L and R registers 
contain the log of the absolute value of X. The log 
is in floating point format. 


Normal Return: 
A normal return is indicated by overflow 'OFF’. 
The common logarithm or (natural logarithm) of the argument 


X is the L and R registers upon return. The log will be in 
floating point format. 


Memory Requirements 
Program: see program listing 

Temporary Space: 13 core memory locations 

Restrictions 

This program uses the following floating point subroutine: 


Name ID Number 
Double Precision Divide 2186650 





Timing 
Approximate Timing: 


FLN - 32.6 milliseconds 
FLG - 36.0 milliseconds 
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Floating Point Common Logarithm or (Natural Logarithm) 
Card Deck 





Symbolic source card deck for this program is 2186647-1 


Appendices 


Flow Chart 


Flow chart for this program is in Appendix A. 


Tables 


Program Listing 


A listing of this program showing memory locations (relative 
to zero), octal representation and source card contents (in 
symbolic language format) for each memory location used by 
the program is in Table 1. 
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Ident: 





ID: FEX 
Introduction 


The purpose of the Floating Point Exponential, E*, program 
is to compute the floating point exponential E* of a 
floating point argument, X. 


The program conforms to the specifications of the Floating 
Point System 2186654. 


This mathematical routine has been prepared to fulfill the 
requirements of Specification Number CDV 86648. This 
drawing, together with the reference documents, describes 
the procedures for preparing program card decks and the 
method of utilizing same. 


The program flow chart shown in Appendix A describes the 
logic used to develop the program. 


Mathematical Method 


To compute E*, the argument X is first reduced to a fixed 
point number. 


x= 2M 
X= M) where M, is a fixed point number. 
Change the exponent of E to an exponent of 2 


p= Alga) 





= 22 +3 where M, is an integer and |M,J <1. 


EX = M2. gltg(1n 2) 
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Mathematical Method (Cont'd) 
Let f= 
Then EX = 22. 


Evaluating Ef, 


fee te andreas et [ay 2 
(D + £°) 





2.000000000000575924 
0049962489 136450764 
4.903154798968682648 
42.013532895 041661680 


For positive xX, 


eX = M2. 





For negative X, 


EX = M2. (2) 
E 


Program Operation 
Program Assembly and Loader 


This program is assembled by the SLAP 2 System (2186651) 
using the symbolic card deck for this program and the 
associated subroutine (section 7.0). The SLAP 2 Loader 
uses the object program card deck, a product of the SLAP 2 
assembly, to load the program into core memory. 


Program Control 


This program is used as a subroutine by a control program. 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS | SIZE | CODE IDENT NO. 


ARE IN INCHES AND INGLUDE TMICNNESS 2/8 6643 
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Control Program 


With the floating point argument X in the left (L) and 
right (R) accumulators, 





Operation Address 


a TSP FEX (see below) 
atl Return 


FEX - Symbolic entrance in the Floating Point Exponential, E*, 
subroutine 


Output Format 
Error return: 
An error return is indicated by overflow 'ON'. 


X> 709.98 = (1023) In 2. The L register will contain 
(37777777) upon return. 


Normal return: 

A normal return is indicated by overflow 'OFF'. 

The exponential of the argument X will be in the L and R 
registers upon return. The exponential will be in normalized 
floating point format. 

Memory Requirements 

Program: see program listing 

Temporary Space: None 


Restrictions 


This program uses the following floating point subroutine: 
Name ID Number 
Double Precision Divide DvD 2186650 


UNLESS OTHERWISE SPECIFIED: DIMENSIONS | SIZE | CODE IDENT NO. 


[ARE IN INCHES AND INCLUDE THICKNESS, ]}) / SE Oo OLE 
(OF PLATING. TOLERANCES ARE: 


OK = se 
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Timing 


Approximate timing: 41.4 milliseconds maximum 





Floating Point Exponential, E%, Card Deck 





Symbolic source card deck for this program is 2186648-1. 
Appendices 

Flow Chart 

Flow chart for this program is in Appendix A. 

Tables 

Program Listing 

A listing of this program showing memory locations (relative 
to zero), octal representation and source card content (in 


symbolic language format) for each memory location used by 
the program is in Table 1. 
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Identification 
ID: FIN or (FCT) 
Introduction 


The purpose of the Floating Point Tangent or (Cotangent) 
program is to compute the floating point tangent or 
(cotangent) of a floating point radian argument, X. 


Tie program conforms to the specifications of the 
Floating Point System 2186654. 


This mathematical routine has Leen prepared to fulfill 
the requirements of Specification Number CDV 86649. This 
drawing, togeter with the reference documents, Jescribes 
the procedures for preparing program card decks and the 
method of utilizing same. 


The program flow chart shown in Appendix A describes 
the logic used to develop the program. 





RADIO CORPORATION OF AMERICA 


ARE IN INCHES AND INGLUDE THICKNESS a] 
OF PLATING, TOLERANCES ARE: 
Ka HO = ii 


FRACTIONS =:t  ANGLES=+1/2° 
DEL 1263-3A (2/65) 








3.0 Mathematical Method 


To evaluate tan X, the following procedure is followed: 
x (2) = 21+ K+ 2 where Izl <o 
Kl =0or1 
I is an integer 
for |z|< .5, set S = Oandw= |2| 
IZ] = .5, sete S=landw= (1 - kl ) 
X<0,K=lsetJ=1 
X20, K=0 set J=1 
X<0,K=O0set J=-1 
X20, K=1 set J=-1 
S=1,K=lorS=0,K=0 
X= wW/S (w?) a 
1, K=0orS=0, K=1 
X= S(W)/w.s 


Where 
sW)-atw {b+ 


a = 0.636619772367597 
b 0.07319486991705 
c + 885605722768290 
d 14..870268625197861 
e 57 .818691387368667 
£ 9.32191895364603 


Cotangent is evaluated by the following relationship. 





Cotangent X = 1/tangent X 
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Program Operation 


Program Assembly and Loader 

This program is assem»led ty the SLAP 2 System (2186651) 
using the symbolic card deck for this program and the 
associated su:routine (Section 7.0). The SLAP 2 Loader 


uses the object program card deck, a product of the 
SLAP 2 assembly, to load the program into core memory. 


Program Control 


This program is used as a subroutine by a control program. 


Control Program 


With the floating point radian argument, X, in the left (L) 
and right (R) accumulators, 


Location Operation Address 





TSP FTN or (FCT) 
(see below) 


atl Return 


Symbolic entrance in the Floating Point Tangent 
subroutine 


Symbolic entrance in the Floating Point Cotangent 
subroutine 
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Output Format 


Error Return: 


An error return is indicated by overflow 'ON'. 


a. The [x] > 26 radians. The L register will be 
positive upon return. 


Te iKkl< gt ont: attempting to evaluate cotan, 


an error return is given. The L register will be 
negative upon return. 


If |x| = N. 7/2, where N is odd, while evaluating 
tan X, an error return will be given. The L regis- 
ter will be negative upon return. 


If |x| = N. 17/2, where N is even, while evaluating 
cotan X, an error return will ‘e given. The L regis- 
ter will be negative upon return. 


Normal Return: 
A normal return is indicated »y overflow 'OFF'. 
The tangent or (cotangent) of the argument X will de 
in the L and R registers upon return. The tangent or 
(cotangent) will be in normalized floating point format. 
Memory Requirements 


Program: see program listing 


Temporary Space: 13 core memory locations 





Restrictions 
This program uses the following floating point subroutine: 
Name Number 


Double Precision Divide 2186650 
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Timing 


Approximate timing: 41.2 milliseconds 


Floating Point Tangent or (Cotangent) Card Deck 





Symbolic source card deck for this program is 2186649-1. 
Appendices 

Flow Chart 

Flow chart for this program is in Appendix A. 


Tables 


Program Listing 


A listing of this program showing memory locations (relative 
to zero), octal representation and source card content (in 
symbolic language format) for each memory location used by 
the program is in Table 1. 
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STORE IN C1 
SBA (C-1) 








STORE INC 
AND C2 
RBA (46+1-c) 
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GET MANTISSA| 
M AND SCALE {¢———— 


AT ZERO 


Bs 


™(?) SHIFT KINTO 
7 sicn Bit AND 
stone 


q [PARRY ouT 






































SCALED AT 
ZERO 











[CALCULATIONS 
[ro compuTE 
PROPER VALUE] 

















SCALE W AT | 
L__J 
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ONE 





FPT/C2 
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EXCLUSIVE OR K BIT 
AND SIGN OF 
RESULT NEGATIVE? 
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Identification 
ID: DVD 


Introduction 





The purpose of the Double Precision Divide program is to 
compute the double quotient of the double precision dividend 
R and the double precision divisor S. 


The program conforms to the specifications of the Floating 
Point System 2186654. 


This mathematical routine has been prepared to fulfill the 
requirements of Specification Number CDV 86650. This 
drawing, together with the reference documents, describes 
the procedures for preparing program card decks and the 
method of utilizing same. 


The program flow chart shown in Appendix A describes the 
logic used to develop the program. 


Mathematical Method 


For ease of notation, let M designate the most significant 
part and L the least significant part of the divisor S. 


2 
“e[-@-@] 
Program Operation 
Program Assembly and Loader 
This program is assembled by the SLAP 2 System (2186651) 
using the symbolic card deck for this program. The SLAP 2 


Loader uses the object program card deck, a product of the 
SLAP 2 Assembly, to load the program into core memory. 
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Program Control 


This program is used as a subroutine by a control program. 


Control Program 


With the dividend R in the left (L) and right (R) accumulators, 


Location Operation Address 





TSP DVD (see below) 
5, M,T (see below) 
Return 

DVD - symbolic entrance in the Double Precision Divide 

M is the double precision divisor S 

T (04 T 6) is the index tag modifying M 

Output Format 

Error Return: 


An error return is indicated by overflow 'ON'. 


The dividend is larger than, or equal to the divisor. The 
difference, (S - R), is the L and R registers upon return. 


Normal Return: 


A normal return is indicated by overflow ‘OFF’. 





The quotient, (R +S), is in the L and R registers upoa 
return. The qu>tient is scaled at (J - K) where J is the 
original scaling of the dividend and K is the original 
scaling of the divisor. 
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Memory Requirements 


Program: see program listing 


Temporary space: none 


Restrictions 


The absolute value of the dividend must be smaller than that 
of the divisor. 


Timing 


Approximate Timing: 7.9 milliseconds 


Double Precision Divide Card Deck 





Symbolic source card deck for this program is 2186650-1 
Appendices 

Flow Chart 

Flow chart for this program is in’ Appendix A. 

Tables 

Program Listing 

A listing of this program showing memory locations (relative 
to zero), octal representation and source card content (in 


symbolic language format) for each memory location used by 
the program is in Table 1. 
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STORE SIGN 
ae IS S NEGATIVE? 
i DIVIDEND, R 
Ps ’ 


IS R NEGATIVE? 











ICOMPLEMENT| 








[COMPLEMENT 
R 








NO 








PICKUP ADDRESS: ¥ 
OF DIVISOR. 
BUILD LDB M,T COMPUTE 


COMMAND. YS 
Ms 

















a 





NORMALIZE MS. 
MS SCALED AT -N. 
N IS NUMBER 

OF SHIFTS 














a ee 


STORE SIGN 
OF s 





| NORMALIZE LS. 
18 SCALED AT -M. 
M IS NUMBER 

OF SHIFTS 

















i N.B: 
Let MB denote most 
i significant part of 8; 
18, least significant 
i an 
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COMPUTE 


BUILD SBA COMMAND R (4)-@) 
Is ms) ~ (ms 
TOSCALE 45 SCALED AT 


@2+M-N) (Q-M+N) 











@ 
SHIFT LS RIGHT 
JONE. LS SCALED 
AT (-M-1) 





EXECUTE SBA| 
COMMAND |, 


TO SCALE 52 











i to 


B ‘COMPUTE 
COMPUTE 5 mR si (28) 215 
SCALED AT ‘3 Ms’ MS |\Ms/ ~ M5] 
1-M+N SCALED AT 0 








¥ 
COMPUTE 


SCALED AT 
2(1-M+N) 





‘SIGN OF DIVIDEND = 
SIGN OF DIVISOR? 











¥ ICOMPLEMENT| 
EXECUTE R 
SBA COMMAND 


5 
1s 
roscate tS 
‘ a) 
COMPUTE 
1s)? _ 8) 
(i) ~ i 


SCALED AT 
Q-M+ny 
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